Possibilities for the synthesis of superheavy element Z = 121 in fusion reactions* 
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Based on the dinuclear system model, the calculated evaporation residue cross sections match well with the 
current experimental results. we systematically study the synthesis of superheavy element Z = 121 through 
combinations of stable projectiles with Z = 21—30 and targets with half-lives exceeding 50 days. The influence 
of mass asymmetry and the isotopic dependence of projectile and target nuclei are investigated in detail. The 
reactions 7°*Es (*°Ti, 3n) 27°7121 and 7°?Es (*°Ti, 3n) 7°°121 are found to be experimentally feasible for 
synthesizing superheavy element Z = 121, with the maximal evaporation residue cross sections of 6.619 fb at 


219.9 MeV and 4.123 fb at 223.9 MeV. 
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I. INTRODUCTION 


The production of new superheavy nuclei (SHN) represents 
a challenging frontier in the realm of low-energy nuclear re- 
action. Over the years, experimental and theoretical nuclear 
physicists have explored SHN synthesis since the prediction 
of the “island of stability” around Z = 114, N = 184 [1, 2], 
and the Skyrme-Hartree-Fock method considers Z = 120, 124 
or 126 and N = 172 or 184 as magic numbers [3]. The syn- 
thesis of superheavy elements (SHEs) Z = 107-112 was ac- 
complished in GSI using cold fusion reactions with Pb and Bi 
targets [4]. However, despite the successful synthesis of SHE 
Z = 113 via cold fusion reaction Zn + ?°9Bi at RIKEN [5], 
the evaporation residue cross section (ERCS) cgp was only 
0.03 pb, reaching the limit of experimental detection at that 
time [6]. To overcome this challenge, researchers in Dubna 
turned to hot fusion reactions with **Ca beams and actinide 
targets. This method results in the production of SHEs Z = 
114-118 [7-12], which completes the seventh period of the 
periodic table. 

In recent years, many new isotopes with Z < 118 were 
synthesized using modern accelerators, such as DC-280 and 
U-400 of Dubna SHE factory, RILAC of RIKEN, SFC of 
HIRFL and UNILAC of GSI [7, 12—16], but the production 
of SHEs Z > 118 remains a challenge. Previous attempts 
to produce SHE Z = 120 using °°Fe+?4*Pu [6] and °4Cr + 
248Cm [17] reactions at Dubna and GSI, respectively, did not 
observe any a decay chains associated with this element. The 
three events reported by GSI in Ref. [17] were later found 
to be random events [18]. In 2020, with the gas-filled re- 
coil separator TASCA at GSI, search for synthesizing the 
SHEs Z = 119 and Z = 120 was conducted via the re- 
actions °°Ti + 24°Bk and °°Ti + 24°CE£, yet neither was de- 
tected [19]. In 2022, RIKEN estimated the optimal incident 
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energy for synthesizing SHE Z = 119 through the reaction 
51v4?48Cm [20]. Therefore, the synthesis of SHEs Z > 118 
requires not only more advanced detection and identification 
techniques, but also an appropriate reaction system. 

In order to describe the process of fusion-evaporation reac- 
tions accurately, several models and different fusion mecha- 
nisms were proposed. The improved quantum molecular dy- 
namics (ImQMD) model [21], time-dependent Hartree-Fock 
theory [22-25], fusion-by-diffusion model [26], cluster dy- 
namical decay model [27], two-step model [28, 29], dinuclear 
system (DNS) model [30—44] and other methods [45—48] are 
proved to be reliable in reproducing experimental data and 
also have given predictions about the synthesis of unknown 
nuclei [22, 45, 49-56]. 

The synthesis and decay of elements Z = 119 and Z = 
120 has been extensively studied theoretically [22, 33, 35, 
45, 52, 57-59], whereas only a limited number of calcula- 
tions have been conducted on the synthesis of SHE Z = 121. 
To address this research gap, this paper aims to investigate 
the optimal projectile-target combinations to synthesize SHE 
Z = 121 and provide a reference for future experimental at- 
tempts. 

The paper is organized as follows. In sect. II, we describe 
the DNS model and examine the predicting reliability. The 
ERCSs of Z = 121 isotopes in different reaction channels 
are discussed in sect. III. Finally, we make a conclusion in 
sect. IV. 


Il. THEORETICAL DESCRIPTIONS 

In the DNS model, the ERCS of synthesizing SHN in the 
center of-mass frame can be obtained with the following ex- 
pression: 


(2J +1)T(Ecm., J) 
7 (1) 


x Pon(Eemm.; J) Weur( Eom. J). 


Here T (Eec.m., J) represents the transmission probability 
for colliding system to overcome the Coulomb barrier Vj. 


e Pon (Ec.m., J) is the fusion probability for forming a com- 
s pound nucleus [60]. Weur (Ec.m., J) denotes the survival 
70 probability that the excited compound nucleus emits neu- 
7 tron, instead of undergoing fission, to reach the ground state 
72 [61]. The nucleus-nucleus interaction potential considering 
73 the quadrupole deformation is expressed as following [62]: 


= 5011 — 68)? + TOC — 68)? 
+ Ve(R, 21, B2, 01,02) 
+ W(R, b1, P2, 1, 62). 


V(R, 81, B2, 01,02) 


(2) 
7s In this formula, 812 and 6f, denote the dynamical 
7s quadrupole deformation parameters and static deformation 
77 parameters for projectile and target nucleus, respectively. 01,2 
7s are the collision angles for deformed projectile and target nu- 
cleus, respectively. The stiffness parameters C1,2 is written 
as [63]: 
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A = 2 represents the quadrupole deformation. The 
ss Coulomb potential Vc is determined using Wang for- 
s4 mula [64]: 
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(4) 
ss The nuclear potential Vy is given by the Woods-Saxon po- 
s7 tential [64]: 
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(5) 
In the capture process, the transmission probability 
T(Eem., B, J) is described via the Ahmed formula [65, 66]: 


1 — exp(—4za) 
T(Eem., B,J) = : 6 
? (Bc. ) 1+ exp(27(87 — a)) (6) 
ə Here a = Ven Hem: Cn and Jj = 


TOTI 

24 RŽ (J) 
93 i aM. H represents the reduced 
94 mass and ay denotes the Morse parameter [67]. 
ə Taking into account the barrier distribution function f (B), 
æ T (Ec.m., J) is written as: 
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The asymmetric barrier distribution parameters are given 
in Ref. [68]. The capture cross section Cecap is calculated as 
follows [62]: 
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Ocap (EBe.m.) = 5 (2J + 1) T 
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The fusion of the colliding nuclei is determined by the po- 
tential energy surface expressed as [62]: 


U (Ni, 21, N2, Z2, R, b1, b2) 
= Er (N1, Z1) + Ep (No, Z2) — Ep (N3, Z3) 
+ Von (M1, Za, No, Z2, R, b1, B2). 


(9) 


The data of the binding energies of the colliding nu- 
cleus Ep (N12, Z1,2) and the formed compound nucleus 
Es (N3, Z3) is taken from Ref. [69]. Von denotes the 
nucleus-nucleus interaction potential. 

The nucleon transfer is treated as a diffusion process at the 
lowest point of the potential energy surface, known as the 
driving potential [62]. To form a compound nucleus, the din- 
uclear system must surpass the inner fusion barrier Bfus along 
the mass asymmetry degree 7 = (Ap — Ar) / (Ap + Ar), 
which denotes the potential energy disparity between the in- 
cident point and the Businaro-Gallone (B.G.) point (the peak 
of the driving potential) [70], defined as Brus = U (nB.c.) — 
U (m). The fusion probability Pon(Ecm., J) is determined 
through the summation of the distribution probabilities of 
crossing the B.G. point P (Ni, Z1, £1, t) as follows: 


Nec. Zea. 
=> 5 PI (Ni, Zi, E1,t = Tim(J)). 


Nı=1 Z,=1 
(10) 
Here the interaction time 7int(J) is calculated via the de- 
flection function method [71]. P (Ni, Z1, Æ1, t) is calculated 
through solving the two-dimensional master equation: 
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adP(M1, 21, F1,t) 


dt 
= >». Wy, ,21:N1,Z1 (t) 


Ni 

x [dn,,z,P(Ni, Z1, E1, t) — 
+X Wm, Zun, (t) 

Zi 

x [dn,,z,P(N1, Zi, E1,t) — dn, z; P(M1, Za, E1,t)] 
— [Aqr(O(t)) + Ass (OJP, Z1, E1, t). 


dn: z, P(N1, Zi, E1,t)] 


(11) 

Here Wy, ,z,; NiZy denotes the mean transition probability 
from state (N1, Z1) to state (Nj, Z1) [72], dy,,z, is the mi- 
croscopic dimension of state (N1, Z1). The quasi-fission rate 


ize Age and fission rate Ags can be given by the one-dimensional 
122 Kramers formula [73]. 

130 The survival process is primarily determined by the ri- 
131 valry between fission and neutron emission [74]. The survival 
132 probability at the excitation energy Ečy can be expressed as: 


Weurl( Een, T, J) 
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: P(E}, J) 
= P(Ety, 2, J) || ees 


i=1 k 


a P (Ečy, 2, J) denotes the realization probability for emit- 
ting x neutrons [75]. E* represents the excitation energy of 
the compound nucleus which have emitted 2— 1 neutrons [57]. 

The neutron decay width r,a (ÆEž¥, J) is calculated using the 


Weisskopf-Ewing theory [76]: 
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Here J, = [0, EX — Ba- ô — +}. 6 and B, represent the 
141 pairing correction and the the neutron separation energy [33], 
142 respectively. The level density p is expressed as in Ref. [77] 
143 and Giny denotes the inverse reaction cross section [78]. 
T;(E3,J) is the fission decay width given by Bohr- 
145 Wheeler transition-state method [79]: 
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vz 1.1A/12 [80, 81]. The temperature-dependent fission barrier 
19 Be(E¥, J) is calculated by the following expression [82, 83]: 


E* 
M * O eas A 
+ BP (E; =0,J)exp ( =) (15) 


(55 ) su + 1). 


Here Be” denotes the macroscopic part of the fission bar- 
ise rier. T; and zup represent the nuclear temperature and the 
1s3 temperature dependent parameter, respectively [82]. Be 
154 is the microscopic shell correction energy in the ground 
155 State [69] and Ep = 25 MeV [50]. Jg.s. and Js.q, are ex- 
156 pressed as in Ref. [84, 85]. 

To evaluate the accuracy of employing our model for pre- 
158 dicting the ERCSs of SHN, we present in Fig. 1 the compar- 
159 isons between calculated ERCSs and the experimental data 
160 in the reactions “8Ca + 24°Cm [12, 86], 48Ca + 248Cm [87], 
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48Ca + 249Bk [88] and 45Ca + ?49Cf [12, 89, 90]. The calcu- 
lation uncertainties arise from the relatively subjective choice 
of Ep range [91]. The fission barrier heavily relies on the 
contribution of the shell correction energy, and the reduction 
of the shell correction energy with increasing excitation en- 
ergy is described by the Æp values, which lies in the range 
10 MeV < Ep < 30 MeV [92]. 

As shown in Figs. 1(a)-(d), the ERCSs show a decreasing 
trend with increasing proton number of compound nucleus. 
For the reactions *8Ca + 24°Cm and *8Ca + ”49Cf, the maxi- 
mal ERCSs of both calculation and experiment appear in the 
3n-emission channels. The 4n-emission channels are more 
favorable for the synthesis of SHN with the reactions “8Ca + 
248Cm and 48Ca + 74°Bk. The predicted ERCSs align well 
with the experimental results, especially for the reaction **Ca 
+ 249Cf, A maximal ERCS of 0.42+9:$7 pb for the reaction 
48Ca + °49Cf is predicted at the 3n-emission channel at Ei, 
= 32.0 MeV, which is consistent with the experimental value 
of 0.5755 pb with E&éy = 32.1 — 36.6 MeV at the same 
channel. This validates the applicability of the DNS model 
in the prediction of synthesizing new elements via fusion- 
evaporation reactions. 


Ill. RESULTS AND DISCUSSION 


To avoid facility contamination by unstable beam, we 
choose stable projectiles with Z = 21 — 30 and actinide 
targets with half-lives exceeding 50 days for the experimen- 
tal duration, and the optimal reaction systems are summa- 
rized in Table 1. One can see that the most favorable re- 
actions and ERCSs (the optimal F..m.) for producing iso- 
topes 29°—302191 are 252Es (46Ti, 3n) 299121, 4.123 fb (223.9 
MeV), 248Cf£ O°V, 3n) 29°121, 0.566 fb (239.1 MeV), 2°4Es 
(46Ti, 3n) 297121, 6.619 fb (219.9 MeV), 2°4Es (4’Ti, 3n) 
298121, 1.331 fb (222.3 MeV), Pm (Sc, 3n) 299121, 
8.778 fb (213.6 MeV), 7°4Es (4°Ti, 3n) 300121, 0.453 fb 
(228.5 MeV), ?54Cf (V, 3n) 391121, 3.705 fb (229.0 MeV) 
and ?54Cf (51V, 3n) 302121, 0.524 fb (234.1 MeV). 

As mentioned in the previous paragraph, the largest maxi- 
mal ERCS corresponding to the synthesis of SHE with Z = 
121 is 8.778 fb in the reaction “°Sc+?°’Fm. Besides, the reac- 
tions *°Ti+?52Es and *°Ti+?°4Es offer large maximal ERCSs 
of 4.123 and 6.619 fb. Considering the experimental fea- 
sibility, The ?54Es target is currently available among sev- 
eral Es targets in laboratory [94], with a half-life of 275.70 
d. The ?°7Es target has a comparatively longer half-life of 
1.29 y, making it a potential target for experimental purposes. 
Therefore, despite a slightly higher ERCS of the reaction 
45Sc+?57Fm, the reactions *°Ti+?°:254Es are more feasible 
for experimental purposes. 

In Figs. 2(a)-(c), we present the calculated ERCSs of the re- 
actions “°Sc+257Fm, 48Ti+254Es and ®!V+2°!Cf. These reac- 
tions yield the same compound nuclei of 9°7121. Notably, our 
analysis reveals a consistent diminishing trend in the maximal 
ERCSs for synthesizing the same isotopes, 799121 via the 3n- 
emission channel and 79°121 via the 4n-emission channel, as 
the charge number of the projectiles increases. This trend can 
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Fig. 1. (Color online) Comparison of the predicted ERCSs with the experimental results [12, 86-90] for the synthesis of Lv (a, b), Ts (c), 
and Og (d). The calculated ERCSs in the 2n, 3n, 4n and Sn-emission channels are denoted by the dashed, solid, dash-dot and dotted lines, 
respectively. The shades indicate the uncertainties of calculated ERCSs. The experimental results for the 2n, 3n, 4n, and Sn-emission channels 
are denoted by inverted triangles, circles, squares and triangles, respectively. 
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Fig. 2. (Color online) The predicted ERCSs of the reactions *°Sc+?°" Fm, **Ti+?°4 Es and °'V +79! Cf. The 3n and 4n-emission channels 
are indicated by the blue solid and red dotted lines, respectively. The shades indicate the uncertainties of calculated ERCSs. 


be attributed to the reduced fusion probability resulting from 
the increased mass asymmetry. To further investigate the in- 
fluence of mass asymmetry on the fusion-evaporation reac- 
tion, the fusion probabilities and the driving potential for the 
reactions “°Sc+?5"Fm, *®Ti+?54Es and °!V+?5!Cf are pre- 
sented in Fig. 3 and Fig. 4. 


Fig. 3 reveals an ascending trend of the fusion probabili- 
ties with the increasing EKé.. This occurs due to the height- 
ened dissipated energy within the dinuclear system at higher 


225 Ečy, thus rendering the reaction system more likely to over- 
22 come the inner fusion barrier. Besides, it can be observed 
27 in Fig. 3 that the reaction *°Sc+?°’Fm exhibits much larger 
228 fusion probability compared to the other two reactions. Con- 
29 versely, the fusion probability of the reaction °!'V+?°!Cf is 
230 the lowest. This significant difference can be attributed to the 


2 


a 


1 different Bfus values influenced by the changed mass asym- 
232 metry. 


23 Fig. 4 reveals that as the mass asymmetry of the reac- 


TABLE 1. The favorable reaction systems for producing SHEs Z = 121. The isotopes, the reaction systems and the half-lives of correspond- 
ing targets [93] are presented in columns 1-3. The optimal incident energies Em, and the EGy are listed in columns 4-5, respectively. The 
maximal calculated ERCSs for certain neutron emission channel are shown in the columns 6. 


Isotope Reaction Ty /2 (target) Ec.m. (MeV) Eén (MeV) CER (fb) 
795191 252 Fs (46Ti 3n) 1.29 yr 223.9 36.0 4.123132. 
296121 2A8CREOV, 3n) 333.50 d 239.1 36.0 0.56646:542 
aero 254F (467) 3n) 275.70 d 219.9 35.0 6.61919 025 
49CFP1V,3n) 351.00 yr 240.3 35.0 OS06 5 is. 
298121 54Rs(47Ti,3n) 275.70 d 222.3 36.0 t3311 33 
299121 25TEm(f5Sc,3n) 100.50 d 213.6 36.0 S77. a. 
254 Fs (48 Ti,3n) 275.70 d 227.6 36.0 Lori. 
52CfP°V,3n) 2.64 yr 232.3 34.0 L368", 
51CfP1V,3n) 898.00 yr 238.2 35.0 C5405) 2 
300191 254 Fs (49T1,3n) 275.70 d 228.5 36.0 0.45370 395 
| aPC VS) 60.50 d 229.0 33.0 S105" 545 
254Es(50Ti,3n) 275.70 d 232.5 35.0 Ona ae. 
2021321 Che Voi) 60.50 d 234.1 34.0 0:5240:$38 
10° — —— 1 r —— 234 tion system decreases, the entrance channel approaches closer 
_ A: 235 to the B.G. point, resulting in a corresponding decrease in 
10 F 3 2 the Brg value. For the reaction *°Sc+?57Fm, the Brg is 
5 ee oer 237 13.1 MeV, which is lower than the reactions *®Ti+?5+Es 
10 F a = F 238 (Bfus=17.1 MeV) and °!V+25!Cf (Brus=17.8 MeV). Conse- 
z 10%} fe ] 239 quently, the reaction 4°Sc+?°’Fm is more likely to overcome 
£ E ef 459 ¢4257F 24 the inner fusion barrier, resulting in an enhanced fusion prob- 
107} if aa aoe 21 ability in Fig. 3. Evidently, the heightened fusion probabili- 
| * 48Tj+Es | 22 ties, stemming from the reduced mass asymmetry, establishes 
a i a 51471251 243 the superiority of Sc and Ti-induced reactions for producin 
10 F ; VRIICE ] p y p g 
a ea SHE with Z = 121. 
10 
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Fig. 3. (Color online) The calculated fusion probabilities of the re- 
actions *°Sc+?°"Fm (black solid line), “*Ti+?°*Es (red dashed line) 
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and °!v+?>!C¢ (blue dash-dot line). 
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(Color online) The driving potential for the reaction 
45 304257 Rm, 1STi+?"*Es and 5tV+°51Cf as a function of mass 2 


asymmetry. The arrows indicate the entrance channel. 
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In Fig. 5, we present an analysis of the calculated maxi- 
mal ERCSs, the corresponding incident energies and Q values 
on the reactions involving 4°—°°Ti projectiles and 7°2:2°4Es 
targets. Fig. 5(a) reveals that the reactions employing the 
neutron-rich ?°*Es target consistently yield larger maximal 
ERCSs compared to those employing the 7°?Es target. More- 
over, the maximal ERCSs decrease as the neutron number of 
the projectiles increases. Notably, the odd-even effects also 
have an impact on the maximal ERCSs, with even-A Ti pro- 


254 jectiles resulting in relatively enhanced ERCSs. Fig. 5(b) il- 
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lustrates that the optimal incident energies for reactions with 
252Fs target are about 3-4 MeV higher than those with ?54Es 
target. Additionally, the optimal incident energies exhibit a 
discernible increase with evident odd-even effects as the neu- 
tron number of the projectiles increases. 


For all the reactions 46-50Ti+?52Es and 46—5°Tj+254Es, the 
corresponding Eġy fall within the range of 35-37 MeV. This 
range has a limited impact on the optimal incident energies. 
The increasing trend in optimal incident energies can be at- 
tributed to the differences in the Q values. Fig. 5(c) reveals 
that a high neutron excess of the target nuclei enhances the Q 
values of the reaction system, while a high neutron excess of 


67 the projectile nuclei exerts the opposite influence. The odd- 


even effects of the projectiles also have a significant influence 
on the Q values, with reactions utilizing the even-A Ti pro- 


270 jectiles displaying relatively suppressed Q values compared 
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Fig. 5. (Color online) (a) The calculated maximal ERCSs, (b) the 
corresponding optimal incident energies and (c) the Q values of the 
reactions *6~°°Tj+?5?Es and 46—°0T}+7>4Es, 


to those employing the odd-A Ti projectiles. 


To investigate the isotopic dependence on the maximal 
ERCSs and corresponding optimal incident energies, a com- 
prehensive investigation of the capture, fusion, and survival 
stages is essential. In Fig. 6(a), we present the capture cross 
sections for the combinations involving *°~°°Ti projectiles 
colliding with ?5?Es and ?54Es targets at excitation energy 
of EGy = 35 MeV and 50 MeV. Notably, the capture cross 
sections exhibit a rising trend with the increase of Eén:; as 
the ability of surpassing the Coulomb barrier increases with 
higher Gy. Furthermore, it is evident that the capture cross 
sections of the reactions involving 7°7Es targets are notably 
enhanced in comparison to those with ?°Es targets. Addi- 
tionally, there is an upward trend in the capture cross sections 
with a higher neutron excess of the projectiles. These trend 
can be attributed to the decrease in the Coulomb barrier. 


In Fig. 6(b) the excitation energies associated with the cor- 
responding Coulomb barriers Vp + Q of the corresponding 
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Fig. 6. (Color online) (a) The calculated capture cross sections of the 
reactions “°~°°Ti+?°*Es and “6~°°Ti+?>Es with Eén = 35 MeV 
and Ečy = 50 MeV. (b) The excitation energies of the correspond- 
ing Coulomb barriers Vp + Q of the reactions 46—50Tj+252Es and 
46—50Tj+254Es. 


reactions are plotted. It can be observed that the V, + Q 
values decreases with the rising neutron excess of the pro- 
jectiles. Moreover, the reaction systems with the ?5?Es target 
exhibit lower Vp + Q values compared to those with 2°4Es 
target. Consequently, the reactions involving ?5?Es as the 
target nuclei, coupled with neutron-rich Ti projectiles, have 
an increased likelihood of overcoming the Coulomb barrier, 
thereby enhancing the corresponding capture cross sections. 
Regarding the fusion process, Fig. 7(a) shows the 
fusion probabilities of the reactions *°—°°Ti+?52Es and 
46—50Ti+?54Es at Hey = 35 MeV and 50 MeV. As the in- 
creasing probabilities of overcoming the inner fusion bar- 
rier, the fusion probabilities are amplified with higher Ečy- 
These fusion probabilities exhibit a decreasing trend with the 
increasing neutron excess of the projectiles. Notably, the 
employment of the neutron-rich ?°*Es target leads to a rel- 
ative enhancement of the fusion probability, which can be 
attributed to the reduced inner fusion barrier. In Fig. 7(b), 
the inner fusion barriers of the corresponding reactions are 
presented. It is evident that the Bfus values rise with the in- 
creasing neutron excess of the projectiles and are higher in 
the reaction systems with the lighter ?5?Es target. This can be 
attributed to the increased mass asymmetry from the projec- 
tiles with higher neutron excess and targets with lower neu- 
tron excess, which subsequently enhances the By,, values and 
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Fig. 7. (Color online) (a) The calculated fusion probabilities of 
the reactions “°~°°Ti+?°?Es and 46-°°Ti+?°4Es with En = 35 
MeV and EGn = 50 MeV. (b) The Brus values of the reactions 
46-5074 4252 Re and 46-507} 4.254 p85 


hinders the fusion process. 


In Fig. 8(a), the survival probabilities of compound nuclei 
in the 3n-emission channel for the reactions 4°—5°Ti+?°?Es 
and *6~°°T}+?>4Es at Ey = 35 MeV and 50 MeV are plot- 
ted. Notably, the survival probabilities exhibit a decreasing 
trend as the EGy increases. This is due to the damped shell 
effect at increased KGy, resulting in diminished compound 
nucleus stability. Additionally, the 7°+Es-based reactions ex- 
hibit relatively higher fusion probabilities with an obvious 
odd-even staggering pattern. This can be ascribed to the in- 
fluence of the Bn and Bẹ values for the corresponding com- 
pound nuclei, as plotted in Fig. 8(b) and Fig. 8(c), respec- 
tively. The compound nuclei formed via the even-A projec- 
tiles are more likely to de-excite through neutron emission, 
due to their relatively higher Br values and lower B, val- 
ues, This behavior results in the odd-even staggering in both 
the survival probabilities and the maximal ERCSs of the Ti- 
induced reactions. Furthermore, the combined effect of the 
B, and Bp values contributes to the generally higher survival 
probabilities in the reactions with the 7°+Es target compared 
to those with 7°7Es target. This dual enhancement in both 
the fusion and survival stages highlights the advantage of em- 
ploying the ?°*Es target for the synthesis of the isotopes with 
Z= 121. 
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Fig. 8. (Color online) (a) The calculated survival probabilities 
of compound nuclei in the 3n-emission channel for the reactions 
16-5074 4252 Es and 46-O T1424 Es with Eon = 35 MeV and Eéy = 
50 MeV. (b) The By values and (c) the B, values of the correspond- 
ing compound nuclei. 


IV. SUMMARY 


The calculated ERCSs using the DNS model are assessed 
with experimental results of the reactions *8Ca + 74°Cm, *8Ca 
+ 48Cm, 48Ca + 249Bk, and “8Ca + 749°C. Our analysis in- 
dicates consistency between the theoretical predictions and 
experimental results. Based on the DNS model, we investi- 
gated the synthesis of SHE Z = 121 employing stable pro- 


345 jectiles with Z = 21 — 30 and actinide targets with half-lives 
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longer than 50 days, revealing that this element is expected to 
be produced via the reactions #°Sc+?°’Fm, “Ti + ?°4Es and 
467; + ?5?Es. Considering experimental feasibility, the reac- 
tions *°Ti + 7°4Es and 4°Ti + ?5?Es are more favorable with 
the maximal ERCSs and optimal incident energies of 6.619 
fb at 219.9 MeV and 4.123 fb at 223.9 MeV. 


We investigate the mass asymmetry effect, revealing the 


353 enhanced fusion probabilities for the Sc and Ti-induced re- 


354 actions. 


Additionally, the influence of the Q values, the 


355 Coulomb barriers, the inner fusion barriers, the fission bar- 
356 riers and the neutron separation energies on the isotopic de- 
357 pendence of reactions with Ti projectiles and Es targets are 
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